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The power and efficiency performance analyses and optimization of a Stirling engine with heat resis¬ 
tance, heat leakage, regeneration loss and mechanical losses are carried out in this paper by using the 
combination of Senft’s mechanical efficiency model with finite time thermodynamics analysis method. 
The analytical formulae for indicated power, shaft power, thermal efficiency and brake thermal efficiency 
for the Stirling engine are derived by assuming that the heat transfer at finite temperature difference 
between the heat reservoirs and the working fluid obeys the linear phenomenological heat transfer law. 
The optimal operating regions for shaft power and break thermal efficiency are obtained and the in¬ 
fluences of heat leakage, regeneration loss and mechanical friction losses on cycle performance are 
investigated through numerical calculations. The obtained results are expected to provide theoretical 
guidelines for the optimal design and operation of practical Stirling engines. 

© 2014 Energy Institute. Published by Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine is a simple type of external-combustion engine with good development prospects. It has the advantage of high ef¬ 
ficiency, low vibration levels, simple structure and can run on any combustible fuels [1-3]. So, more and more attention has been paid to 
research work of Stirling engine. Kolin [4] had shown that with perfect regeneration, the thermal efficiency of an ideal Stirling engine can 
attain the maximum, i.e. the Carnot efficiency. 

In recent years, the thermodynamic researches on the performance of the Stirling engine include many aspects. For example, one is the 
mechanical efficiency model developed by Senft [5-12]. Senft 5—12] studied the mechanical efficiency model of heat engine and applied to 
the performance analysis of the Stirling engine. He studied the characteristics of a fair-comparison class of reciprocating engines including 
Stirling engine and found that the ideal Stirling engine has the maximal mechanical efficiency. Thus taken together, the Stirling engine is 
seen to have the maximal brake-thermal efficiency in all reciprocating engines. Another is the finite time thermodynamic theory [13-22] 
developed by Andresen et al. Meanwhile, much work has been carried out on the finite time thermodynamic performance of the Stirling 
engine cycle and many new findings have been obtained [23-41 . Endoreversible Stirling engine model with heat resistance loss and 
irreversible Stirling engine model with losses of heat resistance, imperfect regeneration, heat leakage and internal irreversibility were built. 
The power, efficiency, exergoeconomic performance, ecological criteria and entropy generation of the common or solar-driven Stirling 
engines coupled to infinite or finite heat-capacity reservoirs with general working medium and quantum working medium were optimized. 

Senft [39] examined the power and efficiency performance of the Stirling engine with heat resistance, internal heat leak and mechanical 
losses, in which heat transfer obeys the Newtonian heat transfer law q°c A(T) by using the combination of his mechanical efficiency model 
[6-12] with the irreversible model in finite time thermodynamics [25,27,29,34 . The irreversibility of heat transfer is one of the most 
important facts in finite time thermodynamics. Many researchers [40-45] have shown that heat transfer law has a great influence on the 
optimal performance of endoreversible and irreversible Carnot and Stirling engine cycles. On the basis of Refs. [12,35,39 , a further step 
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List of symbols 

A 

volume compression ratio 



V 

efficiency 

B 

buffer space 

T 

cycle period, s 

C 

specific heat, W K 



Ci 

associated heat leakage coefficient, W K 

Subscripts 

E 

transmission force engine mechanism efficiency 

B 

buffer space 

F 

flywheel 

c 

compression 

Ki 

slope coefficient of temperature-time variation law in 

e 

expansion 


regenerator, K s -1 

F 

flywheel 

M 

mechanism 

H 

hot side/high temperature heat reservoir 

n 

mole number, mol 

L 

cold side/low temperature heat reservoir 

P 

piston/power output of the engine cycle, W 

m 

optimal efficiency at maximum profit rate 

Q 

the amount of heat transfer, J 

max 

maximum value 

R 

regenerator 

me 

mechanical efficiency 

R s 

universal ideal gas constant, J mol -1 I< -1 

opt 

optimal value 

T 

temperature, K 

P 

power/piston 

t 

time duration of the process, s 

P 

pressure 

V 

volume, m 3 

R 

regenerator 

W 

work, J 

RL 

regeneration heat loss 

X 

working fluid temperature ratio 

s 

shaft 

Xm 

the optimal working fluid temperature ratio at 

t 

thermal 


maximum power output and maximum thermal 

0 

ambient 


efficiency 

1 

isothermal expansion process/hot side working fluid 

* 

x m 

the optimal working fluid temperature ratio at 

2 

isothermal compression process/cold side working 


maximum shaft power and maximum brake thermal 


fluid 


efficiency 

3 

regenerative heating process 

a,(3 

associate heat conductances, W I< 

4 

regenerative cooling process 


made in this paper is to analyze the shaft power and brake thermal efficiency of the irreversible Stirling engine with heat resistance, heat 
leakage, regeneration loss and mechanical losses, in which the heat transfer between the heat reservoirs and the working fluid obeys the 
linear phenomenological heat transfer law q°cA(T -1 ) by using the combination of Senft’s mechanical efficiency model with finite time 
thermodynamics. The results obtained can provide theoretical guidelines for the design and optimization of practical Stirling engines. 

2. Stirling engine model 

The Stirling engine model working with two constant-temperature heat reservoirs is shown in Fig. 1. The temperatures of the high 
temperature and low temperature side heat reservoirs are Th and Tl. The working fluid is assumed to be the steady-flowing ideal gas. 
Between the cold- and hot-side working fluids, there exists a regenerator R. Heat resistance, heat leakage and regeneration loss are 
considered here. However, some kinds of irreversibilities, such as intern heat rate losses transferred by fluid and solid conduction, pressure 
losses caused by internal friction and fluid escape are not considered in this model. 

2.1. Linear phenomenological heat transfer law 

Heat transfer between the working fluid and the heat reservoirs obeys the linear phenomenological heat transfer law. So, in case of 
isothermal expansion and compression Stirling engine process, the heat transferred into (Qi) the engine from the hot heat reservoir and out 
of (Q 2 ) the engine to the cold heat reservoir per cycle are, respectively 

Qi = cfry = iiR g T] In A (1) 

02 = 0 (V - Tf 1 ) t 2 = nR g T 2 In X (2) 

where a and j3 are the associated heat conductances between the engine and the reservoirs, T\ and T 2 are the temperatures of the working 
fluid in the expansion and compression process, t\ and t 2 are the time duration of the expansion and compression process, n and R g are the 
mole number and the universal ideal gas constant of the working fluid, and A = V 2 IV\ is the volume compression ratio. Here, Qi and Q 2 are 
supposed in absolute values. 

3. Regeneration heat loss 

Because of internal irreversibility, the efficiency of the regenerator is less than unity and the process is not an ideal regeneration process. 
So the regeneration heat loss is determined by [35,45] 
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Fig. 1. Model of irreversible Stirling engine. 


Qrl = nC v (\- VR )(h-T 2 ) (3) 

where C v is molar constant volume specific heat capacity of the working fluid and t]r is efficiency of the regenerator. 

During the regeneration process, the temperature of the working fluid T has a linear variation law with the time duration 42] ; if K\ is the 
slope coefficient, it results: 

dT/dt = iKj (4) 

where J <4 > 0 and is determined by the stuffing material of the regenerator. In Eq. (4), the symbol “±” represents the regenerative heating (+) 
and cooling (-) processes, respectively. 


4. Linear phenomenological heat leakage 

The heat leakage between the heat reservoirs is [46]: 

Q; = C,(V -Th 1 )! (5) 

where Q is the associated heat leakage coefficient and t is the cycle period. It is due to the external conductance between the hot and cold 
heat reservoirs. 


5. Thermodynamic optimal correlation between power output and thermal efficiency of Stirling engine 


Combining Eqs. (l)-(5), the heats supplied by the hot heat reservoir (Qh) and released to the cold heat reservoir ( Ql ) are: 
Qh = Qi + Qj + Qrl 


( 6 ) 


Ql = Q 2 + Qi + Qrl 

The cycle period is t = ti + t 2 + t 3 + £ 4 . Combining Eqs. (1) and (2) with Eq. (4) gives: 
T = nRgI\ In A | nR g T 2 In A | 2 - T 2 ) 


( 7 ) 


a 


( Tr 1 -^ 1 ) p(t 


Ki 


(8) 


Using the first law of thermodynamics and combining Eqs. (6)-(8) yields the power output (i.e. the indicated power) and the thermal 
efficiency of the Stirling engine: 


Pi = (Qh-Ql)A = (Qi -02) A = [nRg(h-T 2 )\nX\/z 


( 9 ) 


Vt - (Qh - Ql)/Qh = (Qi - Q 2 )/(Qi + Qi + Qrl) 

Combining Eqs. (l)-(5) with Eqs. (8)-(10) gives: 


( 10 ) 


Pi = riabaRg In XI<[ /(2 aba + ncR g In XK\ T 2 ) 


( 11 ) 
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nabaRg In XI<^T H T L 

Tit -2abaC[M - ncR g In A/< r 1 r 2 QAr + nabI<^T H T L [E g In X + aC v ( 1 - rj R )] 

where: 


( 12 ) 


AT = T h -T l , a = 1 - x, b = (T 2 - T L )(xT H - 7 2 ), c = T 2 (t h - 6 2 xT L ) + 7 H 7 t (<5 2 x 2 - l) (13) 

and <5 = (a//?) 0,5 , x = T 2 /T 1 is the working fluid temperature ratio. 

In Eqs. (11) and (12), for fixed Th, Tl, a and 0, maximizing P and 7j t respect to T 2 by setting dP/dI 2 = 0 and d^ t /dr 2 = 0 yields the same 
optimal temperature r 2 o pt 


T 2 ,o P t = T H T L [T H (-l + 5x) - 715(1 + 5)}/(t 2 h - 5T 2 ) 

And the corresponding maximum power output P; and maximum thermal efficiency rj t are, respectively: 


(14) 


-nmvRg In XI<^ 

2 (2vT H T L - r + - nmuR In XI<^ T H T L 


(15) 


rj t = naahRg In XI<^T H T L j |2aaQAT(xd 2 + gT H T L + gdj^j - ngR g In XK^T H T L CiAT(gl + dk ) - nahI<^T H T L [P g In X + aC v ( 1 
where d, g, hj , k, l , m, r, u and v are expressed respectively, as: 



(16) 


d = r’7 2 - 57 2 , g = 7 H (r 1 +x) - 7i(x5 + 1), j = T L +xa 2 5~ 2 T H , h = ( 5T L - 7 H )(<T 1 7 H - T L )(xT H - T L ) 2 , 
k = 1 - <5 2 x 2 , l = T h - x 5 2 T l , m = 7^5 _1 - 27 H 7,. + o7 L 2 , r = x7^ _1 - 57 t 3 , v = xT H - T L , u = 5x 2 + 2x + <5 _1 


6. Analysis of mechanical loss [5-12] 

A general reciprocating heat engine can be conceptually represented as Fig. 2 [5,12 . In the diagram, M represents the mechanism, F the 
flywheel and B the buffer space where the inside pressure of the working fluid action on the non-working space side of the power piston pp. 
The source of buffer pressure is due to a special enclosure called buffer space. Usually, the buffer space is open to the surrounding atmo¬ 
sphere and thus the buffer pressure is the atmospheric pressure. The buffer gas absorbs stores and returns energy to the working fluid. It acts 
directly on the power piston and so it diverts and recycles some work W B + away from the mechanism; the flywheel have a similar energy 
interaction W F +. The flywheel and buffer pressure energy interaction reduces the mechanical friction losses of the engine. 

The arrow in the diagram shows how the work transfer between the components. W s is the useful work coming out of engine per cycle at 
the shaft engine. It is the difference between the cyclic work received by the buffer space ( W B+ ) and flywheel ( W F+ ) through the mechanism. 
And the cyclic work taken from flywheel (W F ~) and buffer space (W B ~) through the mechanism sustain the operation of the engine cycle. So, 
the amount of W s depends on the transmission force engine mechanism efficiency (E). E is defined as the ratio of the actual output force to 
the ideal output force for any given input force applied [8,12 . 

The Stirling engine mechanical efficiency is the ratio of this cyclic shaft work W s to output cyclic work W } which is the difference between 
the absolute expansion work W e and the absolute compression work W c of the engine cycle. For fixed volume compression ratio A = V 2 /Vi 
and working fluid temperature ratio x = T 2 /T, the cyclic mechanical efficiency can not exceed to 5—12] : 

Vme = (E,x, A) = E - (E- 1 - E)S(x, A) (18) 



Fig. 2. Conceptual diagram of a reciprocating heat engine [5]. 
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0 xA < 1 

S(x, A) = { x In x - (1 + x) [ln(l + x) - ln(l + A)] - In A 


(1 -x)ln A 


xA > 1 


(19) 


The upper bound ri me {E,x,X) is the mechanical efficiency of an ideal Stirling engine buffered at its optimal constant pressure and having a 
mechanism of constant transmission force mechanism efficiency E. 


7. Shaft power and brake thermal efficiency of Stirling engine 

Combining Eq. (15) of Stirling cycle power output with Eq. (18) of the engine mechanical efficiency gives the shaft power of the engine: 

Ps = PiVme (20) 

When the values of a , 0, T Ht T L , K\, A and E are set, shaft power P s is the function of working fluid temperature ratio x. Setting dP s /dx = 0 
gives the optimal temperature ratio x* m and the corresponding maximum shaft power. In the following, numerical examples are provided to 
analyze the optimal performance of the Stirling engine. In the calculations, ideal gas is chosen as the working fluid. It is set that n = 1.0 mol, 
a = ft = 5.0 x 10 s W K, C v = 20.77 J/(mol K), T H = 800 K. T L = 300 K. K x = 8.0 x 10 3 K/s and X = 4. 

Fig. 3 shows the shaft power P s versus the working fluid temperature ratio x characteristic for different values of E. In the calculations, 
Q = 2.0 x 10 4 W K, t]r = 0.8 and E = 0.8 and 1.0 are set. The curves of P s versus x are parabolic-like ones. From Eqs. (18) and (20), it can be 
inferred that when E = 1.0, the mechanical efficiency ri me = 1.0 and shaft power P s equals the maximum indicated power P/, as shown by 
curve 1 in Fig. 3. It follows from the figure that the maximum shaft power point can never exceed the maximum indicated power point, i.e. 
Ps,max < Pi, max and x* m < x m , with equality holding if and only if 1/A > x m . So, one can be see how mechanical losses reduce the output shaft 
power and decrease the value of optimal working fluid temperature ratio corresponding to the maximum shaft power. 

Combining Eq. (16) of the Stirling cycle thermal efficiency with Eq. 18) of the engine mechanical efficiency gives the break thermal 
efficiency of the engine: 


Vs — VtVme (21) 

Fig. 4 shows brake thermal efficiency r\ s versus the working fluid temperature ratio x characteristic for different values of E. In the 
calculations, it is set that Q = 2.0 x 10 4 W K, tjr = 0.8, E = 0.8 and 1.0. The curves of % versus x are parabolic-like ones. When E = 1.0, the 
mechanical efficiency ^ me = 1.0 and brake thermal efficiency % equals the thermal efficiency rj tt as shown by curve 1 in Fig. 4. It can be seen 
that the maximum brake thermal efficiency point can never exceed the maximum thermal efficiency point. As E decreases, the optimum 
temperature ratio x m corresponding to the maximum thermal efficiency ^t >ma x is slightly higher than x* m corresponding to the maximum 
brake thermal efficiency % >max - One can see that the temperature ratio x m corresponding to the maximum brake thermal efficiency in Fig. 4 is 
different from the temperature ratio x m corresponding to the maximum shaft power in Fig. 3. In Eqs. (11) and (12), T 2 and x are set as two 
variables. Their exist the same optimal temperature T 2 i0p t which lead to the maximum output power P/ and the maximum thermal efficiency 
rj t . Both Pi and 7j t are functions of temperature ratio x. For a heat engine, the maximum power operation regime and the maximum efficiency 
operation regime are usually attained at two different working points, which lead to the different values of x m in Figs. 3 and 4. 

Fig. 5 shows the influence of heat leakage on shaft power P s versus the brake thermal efficiency % characteristic. In the calculations, it is 
set that t]r = 0.8, E = 0.8 and Q = 0,2.0 x 10 4 W K and 5.0 x 10 4 W K. When there is no heat leakage in the cycle (Q = 0), the curve of P s -r]s is a 
parabolic-like one, as shown by curve 1 in Fig. 5. There exists the optimal brake thermal efficiency which leads to the maximum shaft power 
corresponding to the maximum power operation regime. When there exists heat leakage in the cycle (Q = 2.0 x 10 4 W I< and 
Q = 5.0 x 10 4 W K), the curves of P s -ris are loop-shaped ones, as shown by curves 2 and 3 in Fig. 5. Both the brake thermal efficiency and the 
shaft power of the engine have a maxima corresponding to the maximum power operation regime and the maximum efficiency operation 
regime [47 . The maximum shaft power P s>max will not be affected by the heat leakage. 77 s>max and 77 s>opt decrease while P s>opt increases as heat 
leakage increases. 



Fig. 3. Shaft P ower P s versus working fluid tem P erature ratio x for different E. 
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Fig. 4. Brake thermal efficiency rj s versus working fluid temperature ratio x for different E. 





Fig. 5. Influence of heat leakage on the shaft power P s versus the brake thermal efficiency rj s characteristic. 


The optimal operation regions for the irreversible Stirling engine with heat leakage should be located in the regions where the P s -ris 
curve has a negative slope. Thus, the optimal shaft power and brake thermal efficiency of the engine should be: 

Ps,opt < Ps < Ps,max (22) 

and 

Vs,opt < Vs < Vs,max (23) 

where P s ,max, Ps, opt, Vs, max and Tjs, opt are important performance parameters which determine the lower and upper bounds of the optimal 
shaft power and brake thermal efficiency. 

Fig. 6 shows the influence of the regeneration heat loss on the shaft power P s versus the brake thermal efficiency p s characteristic. In the 
calculations, Q = 2.0 x 10 4 W K, E = 0.8 and the efficiency of the regenerator is set as tjr = 0.9, 0.8 and 0.7, respectively. The corresponding 



Fig. 6. Influence of regeneration heat loss on the shaft power P s versus the brake thermal efficiency 7] s characteristic. 
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curves are shown by curves 1, 2 and 3 in Fig. 6. As regeneration loss increases, the maximum efficiency ^ Sima x and the optimal efficiency ri Siopt 
decrease; while the maximum and optimal shaft power P s ,max and P Si0p t will stay unchanged. 

8. Conclusion 

The power and efficiency performance analyses and optimization of a Stirling engine are carried out in this paper by using the com¬ 
bination of Senft’s mechanical efficiency model with finite time thermodynamics analysis method. The model built in this paper represents 
several typical kinds of losses of an ideal Stirling engine: heat transfer loss, heat leakage, regeneration loss and mechanical loss. Specially, the 
heat transfer of the engine model between the working fluid and the heat reservoirs obeys the linear phenomenological heat transfer law 
q°c A(T _1 ). The analytical expressions of shaft power and brake thermal efficiency are derived. Numerical examples show that heat leakage, 
regeneration heat loss and mechanical losses have great influences on the power and efficiency performance of the Stirling engine. The 
optimal operation regions for power and efficiency have been determined. The results can provide some guidelines for the design and 
optimization of the Stirling engine. 
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